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ABSTRACT 

We have performed a sequence of high resolution hydrodynamic simulations of structure formation 
in a ACDM model to study the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects. Including only 
adiabatic gas physics, we demonstrate that our simulations for the thermal effect are converged down 
to sub-arcminute scales. In this model, the angular power spectrum of CMB anisotropies induced by 
the thermal effect peaks at £ ~ 10'*, and reaches an amplitude just below current observational upper 
limits. Fluctuations due to the kinetic effect are a factor of ~ 30 lower in power and peak at slightly 
smaller angular scales. We identify individual SZ sources and compute their counts as a function of 
source strength and angular size. We present a preliminary investigation of the consequences of an early 
epoch of energy injection which tends to suppress power on small angular scales, while giving rise to 
additional power on large scales from the reheated IGM at high redshift. 

Subject headings: cosmic microwave background - cosmology: theory - galaxies: clusters: general - 
large-scale structure of universe - methods: numerical 



1. INTRODUCTION 

On angular scales below 10', the microwave sky carries 
the imprint of large-scale structure in the low-z universe. 
In particular, cosmic microwave background (CMB) pho- 
tons propagating through the universe are inverse Comp- 
ton or Doppler scattered by hot electrons along their path, 
either in dense structures such as clusters of galaxies or 
more generally in hot gas in the intergalactic medium. In- 
verse Compton scattering conserves the number of pho- 
tons but preferentially increases their energy, leading to 
a spectral distortion whose amplitude is proportional to 
the product of the electron temperature and density (or 
pressure). Doppler scattering induces an intensity fluc- 
tuation with the same spectral shape as the CMB itself. 
These effect s were first described by Sunyaev & Zel'dovich 
(1972, 198C) and are known as the thermal and kinet ic SZ 
effects, respecti vely ( for recent reviews see Rephaeli 1995 
and Birkinshaw 1999). 

The thermal effect is one of the primary sources of sec- 
ondary anisotropies in the CMB on small angular scales. 
The change in the (thermodynamic) temperature of the 
CMB resulting from scattering off non-relativistic elec- 
trons is 
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where x = hv/kTcMB — J^/56.85 GHz is the dimension- 
less frequency, and the second expression is valid in the 
Rayleigh- Jeans limit which we shall assume henceforth. 
The quantity y is known as the Comptonization parame- 
ter and is given by 
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where the integral is performed along the photon path. 
Since Tg ^ Tcmb the integrand is proportional to the 
integrated electron pressure along the line of sight. 



The kinetic SZ effect arises from the motion of ionized 
gas with respect to the rest- frame of the CMB. The result- 
ing temperature fluctuation is AT/T = —h, where 
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measures the magnitude of the effect along the line of sight 
if Vr is the radial peculiar velocity of the gas. The different 
dependence on frequency of the two effects can, in princi- 
ple, be used to disentangle them observationally, though 
the primary CMB anisotropies provide "noise" for the ki- 
netic SZ signal. The imprint of the SZ effects on the CMB 
has been s tudied by a numb er of authors (Ostriker & Vish- 



1986t Persi et al. |1995|, d a Silva et al. gOOOj, Refregier 
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200C , among others) 
exhibit substantial quantitative 



et al. |2000| , Seljak et al. 
oretical predictions stil 
uncertainty. 

In this paper, we report first results of a sequence of high 
resolution hydrodynamic simulations designed to study 
the SZ signal on small angular scales. Our simulations can 
be used to calibrate and test the semi-analytic modelling 
that has been done before, explore the effect of varying 
the physics included, and provide predictions for upcom- 
ing experiments. In Section 2, we describe our simulation 
set and we discuss the techniques we use to compute SZ 
maps. In Section 3, we present our results for the power 
spectra of thermal and kinetic SZ effects, and we show SZ 
source counts as a function of source strength and source 
size. Finally, we summarize and discuss our results in Sec- 
tion 4. 

2. METHOD 

We are interested primarily in the consequences of as- 
sumptions about the physical state of the gas for pre- 
dictions of the SZ signal, so we have focussed on one 
particular cosmological model, chosen from the ACDM 
family. For definiteness, we have adop ted th e "concor- 
dance" model of Ostriker & Steinhardt (1995) which has 



Vl^ = 0.3, = 0.7, Hq = lOO/ikms^^Mpc"^ with 



2 



10^ 



10' 



A 

i 10'' 



1 1 

r 1 keV @ z = 7 - ^ 




: 0.3keV @ z = 7 'l^-^,^ " " - ^ 




rO.l keV @ z = 7 - 


- 


- /v 

: 67A / / 




! "^^'^^O' 134A 




1 





injected. The amount of extra energy was chosen to be 
plausible given the observed luminosity-temperatur e rela- 
tion in clusters of gala xies ( Cavaliere, Mend & Tozzi 
Wu, Fabian & Nulsen 1998) and the expected X-ray 
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Fig. 1. — The evolution of the density weighted temperature with 
redshift. The soHd and dashed Hnes are the simulation results, while 
the dotted line shows our predictions based on the Press-Schechter 
theory. 

h = 0.67, f^B = 0.04, n = I and erg = 0.9 (correspond- 
ing to 5h = 5.02 X 10^^). This model yields a reasonable 
fit to the current suite of cosmological constraints and as 
such provides a good framework for making realistic pre- 
dictions. 

We ran a number of simulations on the Cray T3E at 
the San Diego Supercomputing Center using t he GA DGET 
Tree/SPH code (Springcl, Yoshida & White pOOO| ). Each 
simulation employed 2 x 224"^ particles and was carried 
out on 64 processors. Our basic model had a box size of 
134 /i~^Mpc, so the dark matter particles had masses of 
1.5 X IQio h-^MQ and the SPH particles 2.4 x 10^ H'^Mq. 
The SPH densities were computed from 32 neighbours, so 
our minimum gas resolution was roughly 8 x 10^" /i~^Mq. 
Our base simulation (134A) was evolved from z = 50 
to z = 0, and the gravitational f orce softening was of a 
spline form (e.g. Hernquist & Katz 1989| ), with a Plummer- 
equivalent softening length of 25 /i~^kpc comoving. 

As a test of our numerical resolution we also simulated 
the same cosmological model in a 67h~^ Mpc box (evolved 
from z = 100 to the present; 67A), with a softening length 
half that of the 134A run and eight times better mass 
resolution. In all our runs, the full simulation box was 
output at 333 different times, between redshifts z ~ 19 
and 2 = 0. The simulations reported here followed only 
adiabatic gas physics. Since the SZ effect is dominated by 
massive objects with large cooling times this is a reason- 
able first approximation. 

To investigate the effect of an early epoch of energy 
injection we took the z = 7 output of 134A and added 
0.1 keV, 0.3 keV or 1 keV of energy per (gas) particle. This 
injection scenario was deliberately chosen to be somewhat 
artificial. A more realistic approach might be to add the 
energy only to gas within halos above a certain density, to 
describe energy input from galaxy winds or QSO activity 
for example, or to add it over an extended period of time. 
Unfortunately, the very large parameter space of these pos- 
sibilities cannot be explored easily. We thus decided to ex- 
amine an "extreme" scenario to get a handle on the depen- 
dence of the small angle SZ signal on the amount of energy 



ground from the diffuse intergalactic medium (Pen 1999). 
We added the energy at z = 7, which we also took to be 
the epoch of reionization, so that the energy injection takes 
place "before" the gas contributes significantly to the SZ 
signal. A visual inspection of slices through the simulation 
shows that an energy injection as large as 1 keV literally 
blows gas out of small halos, leaving only the gas in the 
larger collapsed objects. The gas is thus more smoothly 
distributed in addition to being hotter. At z = the 
slope of the mass weighted temperature-mass relation in 
the IkeV simulation is thus changed from a canonical 1.5 
to a considerably steeper 1.7. 

The evolution of the density weighted temperature in 
our simulations is shown in Fig. [|. Note the good conver- 
gence between the large (134A) and the small box (67A) 
at redshifts z < 5. In the 67A run, the adiabatic decay 
of the mean temperature is halted at a somewhat higher 
redshift, as expected: due to the better mass resolution of 
this simulation, nonlinear structures of small mass collaps- 
ing at earlier times can be resolved. In the runs with en- 
ergy injection, the temperature increases discontinuously 
at z = 7, and then declines again with the adiabatic expan- 
sion, until shock heating takes over at z ~ 1 — 3, depending 
on the amount of injected energy. At z = 0, all runs yield 
a mean mass- weighted temperature of ~ 0.3 keV. 

We pro duced maps in a manner similar to da Silva et 
al. ( 2000| ). First, the various outputs from the boxes were 
"stacked" back along the line of sight, with each box ran- 
domly translated and oriented in the direction of either 
the X, y or z axis. Then, a grid of 512^ rays subtending a 
constant angle of 1° from the observer was traced through 
the boxes starting at z = 19. For our choice of an open- 
ing angle of 1°, the full 134ft,~^Mpc box just subtends an 
angle equal to the field of view aX z — 19. For the smaller 
67/i~^Mpc run, the field of view reaches the box size at 
a distance z ~ 2, beyond which we covered it by periodic 
replication of the box. 

We produced maps of the y-parameter, the Doppler b- 
parameter, and the projected gas and dark matter densi- 
ties. For example, the y temperature decrement along any 
ray was calculated by distributing the product of pressure 
and specific volume of the gas particles over the angular 
grid: 

= y^PaWg.ij , (5) 

where 'Wa,ij is the value of the (projected) smoothing ker- 
nel (normalized to unity for the pixels covered) of particle 
a at angular grid position {ij), and Lpj^ is the physical 
area of a pixel at the particle's distance. 

The product of pressure and volume, Pa-, for each SPH 
particle is 

p = (7 - 1)(1 - Fp)mM^a;e , (6) 

where 7 is the ratio of specific heats, Yp = 0.24 the primor- 
dial "^He mass fraction, m the particle mass, u the internal 
energy per unit mass, fi the mean molecular weight and 
Xe the ratio of electron and hydrogen number densities. 
In the runs reported here, we do not track Xe dynami- 
cally. Instead, we set it to 1.158 for temperatures above 
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Fig. 2. — Maps of the thermal (top) and kinetic (bottom) SZ 
effects. The maps are 1° on a side and cover both the same field of 
view, here computed for the 134A simulation. Notice that unlike the 
Comptonization y, the kinetic signal b may be negative; we map the 
negative values onto the lower half of a logarithmic grayscale. Struc- 
tures appearing in white are moving towards the observer, those in 
black away. 

10"* K, corresponding to full ionization, and take the gas at 
lower temperature to be neutral. In reality, the inclusion 
of an ionizing background would render this assumption 
invalid at low redshift, but the small densities of the af- 
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Fig. 3. — The angular power spectrum of the SZ effect from an 
ensemble of 15 maps. The upper solid line shows the primary CMB 
anisotropics for this model. The points indicate the mean and error 
on the mean from our ensemble of maps. Filled circles mark results 
for the 134A simulation, open diamonds are for the smaller 67A-box. 
The dotted lines indicate the field-to-field variance between individ- 
ual 1° maps. The arrows at upper right are 95% CL_u{iDer limits 
from a variety of pyppr iments- STT 7TF, (Church et al. 1997), ATC A 
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Ryle (Jones 1997), and 
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19981), BIMA (HplzaDfel et al. 200C), 
(Partridge et al. 1997). 



fected gas do not lead to any significant contribution to the 
SZ signals. These are thus unaffected by our assumption 
about the ionization state of the gas, provided ionization 
happens late enough. On the other hand, it is crucial to 
assume that the gas is neutral at very high redshift, oth- 
erwise the high physical density of the ambient gas in this 
regime would give rise to a significant contribution to the 
SZ effects. In particular, the finite box size of our simu- 
lations would then lead to a strong signal for the kinetic 
effect on large angular scales, reflecting motion of gas due 
to the largest modes in the box. 

3. RESULTS 

We show typical examples for our thermal and kinetic 
SZ maps in Fig. ^. These maps were produced by coadding 
a large number of partial maps, each giving the contribu- 
tion of one of the boxes that we stacked along the photon's 
path. It is interesting to note that filamentary structure is 
easily detected in these partial maps with their small depth 
in redshift space. However, in the full projection along the 
backward photon path, filaments are largely hidden in the 
high level of background arising from the summation over 
many of these structures. Things may not be quite so 
bad for the kinetic effect, where it seems that a somewhat 
larger degree of filamentary structure manages to survive 
at appreciable signal. Note that sources that are "bright" 
in the thermal effect are not necessarily among the bright- 
est in the kinetic effect, and vice versa. Often, the kinetic 
effect shows neighboring large peaks of opposite sign. This 
may be useful in strategies for identifyin g clus ters and su- 
perclusters (Diaferio, Sunyaev & Nusser 2000| ). 
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Fig. 4. — Published estimates of the thermal SZ power spectrum 
compared to the result of this work. Filled circles show our mea- 
surement, while other symbols give comp utations by the following 
authopi^K (Komatsu & Xila^ama |L99E| ), MB (Molpaji^ Birkin- 
shaw pOOq ), CCL4Ckaoray bood), ZH (7Eajag_Xj Pen pOOol), RKSP 
(Refregier pt al gOOo|), SBFTSeljak et al. pOO(f , and : 
Carlstrom 1990- 
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Fig. .5. — Energy injection and its impact on the angular power 
spectrum of the SZ effect. The upper solid line shows the primary 
CMB anisotropics for this model, and the lower solid lines give the 
secondary anisotropics for thermal and kinetic SZ effects in our 134A 
simulation. Dashed, dotted, and dot-dashed lines show the effect of 
an injection of 0.1, 0.3, and 1.0 keV of energy at z = 7, respectively. 
Results are based on ensembles of 15 maps for each of our numerical 
simulations. 



3.1. Power spectra 

For each map we computed a number of statistics, and 
we in general averaged the results for 15 random lines of 
sight to reduce field-to-ficld variance. Of most interest 
here are the angular power spectrum of the thermal and 
kinetic SZ effects (Fig. |). 

A comparison of our basic 67 h^^ Mpc and 134 h^^ Mpc 
runs shows that for the physics we have included we have 
converged in the mean mass-weighted temperature and 
the angular power spectrum for £ < 20000, albeit some 
small residual resolution effects are visible for the kinetic 
effect. The good level of agreement for the two different 
runs down to very small angular scales builds confidence 
in our numerical technique. While this is to be expected as 
the signal is dominated by structures of mass > 10^'^ Mq, 
which should be well covered by our resolution, it also sug- 
gests that fini te bo x size and the box stacking technique 
(White & Hu 2000| ) are not a source of significant uncer- 
tainty in our results. 

We note that the small box gives slightly more power 
at ^ ^ 20000 for the thermal effect, indicating its better 
resolution, while at lower £ the convergence is excellent. 
This is also reflected in the mean Comptonization, which 
both models give as (y) = 2.6 x 10~^. For the kinetic 
effect, we observe a higher level of additional small-scale 
power for the small box on one hand, and slightly less 
power for all £ below £ ~ 20000 on the other hand. This 
suppression likely originates from a deficit of large-scale 
power in the small box, which reduces the amplitude of 
peculiar velocities. 

Interestingly, the power spectrum of the thermal SZ 
effect is not far from current observational upper limits, 
making it a highly interesting target for upcoming exper- 
iments. The power in fluctuations due to the kinetic SZ 
effect is substantially smaller, but may not be out of reach 



in future experiments, especially when multi-frequency ob- 
servations are employed. 

It is interesting to compare published results for the 
power spectrum of the thermal SZ effect with our mea- 
surement, as we have done in Figure ^. There are order of 
magnitude differences between the various results. Note 
that slight differences in the adopted cosmological mod- 
els can partly account for these discrepancies, but even 
for cosmological parameters that are very close, substan- 
tial uncertainties remain. One group of results in this 
comparison is based on analytical computations, using ei- 
ther ext ensio ns of Press-Schechter theory (Komatsu & Ki- 
tayama p^99t Holderfc Carlstrom [T999| , Cooray |200d| , Mol- 
nar & Birkins haw |2000| ) or non-linear perturbation theory 
(Zhang & Pen 2000). Compared to our result, these com- 
putations suggest more power on large angular scales, and 
their spectrum usually peaks at somewhat larger angular 
scale, corresponding to £ of a few thousand. On the other 



hand, the numerical work by Seljak et al. (2000), which is 
based on simu lation s using a moving-mesh hydrodynami- 
cal code (Pen 1998), suggests a very similar shape for the 
spectrum, but with an amplitude whi ch is a factor 1.5-2 
lower than our result. Refregier et al . (|2000| ) use the same 
simulation technique as Seljak et al. ( |2000p , albeit at lower 
resolution and in combination with a different analysis of 
the simulation outputs. Their spectrum peaks at an an- 
gular scale of £ ~ 2000, roughly at the upper boundary of 
the author's own range of confidence of 200 < £ < 2000 
for their result. 

In general, predictions in the literature for the kinetic 
SZ effect show a similar degree of quantitative uncertainty. 
Note however that the computations by Bruscoli et al . 
(|l99|) are very close to our result, as is Hu & White (E99q) 
if reionization is assumed to happen late. On the other 
hand, the signals computed for inhomogcneous rcioniza- 
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Fig. 6. — Thermal SZ sources identified with SExtractor. A circle 
is drawn for each source, with area equal to the area determined 
by the source detecticm software. The map shown in this example 
is the same as in Fig. H, but here displayed with reversed video for 
graphical clarity. 



tion (Gruzinov & Hu |l99g , Benson et al. 200C , Valageas 
et al. ^000| ) are typically an order of magnitude smaller in 
power at the peak of the spectrum. 

The relatively large scatter in these computations high- 
lights that there is still substantial quantitative uncer- 
tainty in the theoretical predictions for the SZ power spec- 
trum. High-resolution simulations like the ones discussed 
here should help to lead to a more reliable answer. 

An additional complication besides computational un- 
certainties is that the theoretical predictions depend quite 
sensitively on the physics of reionization and preheat- 
ing, especially for scales I ^ 2000, where the SZ effect 
dominates over primary CMB fluctuations. This is seen 
when analyzing the SZ signals of our energy injection runs 
(Fig. ^). Depending on the amount of energy injected, 
fluctuations on small angular scales can be suppressed very 
strongly, while additional power on large scales is seen from 
relatively diffuse, ionized gas. This signal on large angular 
scales is generated in a small redshift range right after the 
energy has been injected. Note that on the largest scales 
we probe, the amplitude of the kinetic SZ effect then be- 
comes dominated by finite box size effects as the signal 
reflects the largest modes that are present in the simula- 
tion. 

Apart from influencing the power spectrum, energy in- 
jection also alters the mean Comptonization of the mod- 
els. It becomes (y) — 6.9 x 10~^ for the 0.1 keV model, 
16.6 X 10"^ for the 0.3 keV run, and reaches 38.0 x 10"^ for 
an injection of 1 keV. The latter is thus strongly excluded 
by the upper limit of 15 x10" ^ from the COBE FIRAS ex- 
periment (Fixsen et al. 1996), and even the 0.3 keV model 
is marginally excluded. This highlights the power of the 
FIRAS constraint. 



3.2. Source counts 

The maps for the thermal SZ effect are dominated by 
discrete sources, as suggested by visual inspection. We 
have used t he sou rce detection software SExtractor (Bertin 
& Arnouts 1996) to count individual sources, and to mea- 
sure their strength. A si milar technique has been employed 
by da Silva et al. ( 200C| ). Figure || shows a typical source 
detection, here obtained for the map shown in Fig |[ To 
mark each source, we have drawn a circle with area equal 
to the area of the ellipse matched by SExtractor to each 
source. We used default settings for the detection algo- 
rithm of the image processing software and let it estimate 
background and 'noise' levels automatically. Note that our 
maps are in principle noise-free, so what is interpreted as 
noise in this procedure effectively arises from limitations 
due to source confusion, which leads to non-detections of 
some of the faintest and smallest sources. In a typical anal- 
ysis with SExtractor, we resolve slightly more than half of 
the total SZ signal into discrete sources. 

We measure the strength of a source as the monochro- 
matic brightness change 



Jn Jn 



(7) 



of the CMB integrated over the solid angle of the source. 
Here B,^ is the Planck spectrum of the primary CMB, and 
f{x) is the spectral function 



fix) = X- 
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with X — hv jkToMS,. In the following, we quote results 
for a frequency of 150 GHz. 

In Fig. 1^, we show the cumulative source counts per 
square degree found in the simulations 134A and 67A. We 
find remarkably good agreement between the two simu- 
lations, which differ by a factor of 8 in mass resolution. 
However, the higher resolution run 67A resolves a slightly 
larger number of faint (small) sources, while the bigger 
box 134A shows a small increase of the counts at the 
bright (large) end of the distribution. Both of these weak 
trends reffect expected effects due to box size and numer- 
ical resolution. Notice that the source counts we find are 



quite close to predictions by De Luca et al. (1995) based 
on Press-Schechter theory, however our simulations don't 
show their exact power-law behavior at source strengths 
above 1 mJy. The dashed line in Fig. ^ gives the counts 
for the simulation with 1 keV energy injection. Interest- 
ingly, the counts for the bright sources are unaffected, but 
sources below a flux of 1 mJy (size of ^30") are strongly 
reduced in number density, and still fainter (smaller) ones 
are practically eliminated. 

Finally, we plot in Fig. ^the differential flux distribution 
as a function of source strength. It is seen that the bulk of 
the total flux is provided by sources around 1 mJy, while 
the contribution by very faint and very bright sources is 
negligible. The good agreement at the faint end between 
the two adiabatic runs suggests that both of these simu- 
lations already resolve all the relevant sources that con- 
tribute significantly to the total Comptonization. Consis- 
tent with the results for the source counts, the energy in- 
jection run shows a strong suppression of faint SZ sources, 
while the brightest sources are nearly unaffected. 
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Fig. 7. — Cumulative counts of thermal SZ sources per square 
degree for models 134A (filled circles) and 67A (diamonds); (a) 
as a function of source strength, and (b) as a function of source 
radius. For each source, the source strength Sv was defined as the 
integrated monochromatic flux decrement, integrated over the solid 
angle of the source. The radius of a source was computed based on 
the area detected by the source extraction software. In both panels, 
the dashed line shows the counts for the run with an injection of 
1 keV of energy at z = 7. 



4. DISCUSSION 

We have performed a sequence of high resolution hydro- 
dynamic simulations of structure formation in a ACDM 
model to investigate the thermal and kinetic Sunyaev- 
Zel'dovich effects. Figs. ^ ||, ^ and || show that including 
only adiabatic gas physics our simulations of the thermal 
effect are converged down to sub-arcminute scales. In par- 
ticular, our results are well converged at the peak of the 
power spectrum, and they predict a mean Comptonization 
of (y) = 2.6 X 10^® for the adopted cosmology. 

The detection of the brightest SZ sources with a density 
of ~ 1 per square degree requires an angular resolution of 
~ 1', and the capability to separate brightness fluctuations 
of size ~ 10 mJy (for 150 GHz) from the primary CMB fluc- 
tuations. Several current experiments are capable of meet- 



FlG. 8. — Differential flux distribution as a function of source 
strength for models 134A (filled circles) and 67A (diamonds). For 
each source, a source strength Sv was defined as the integrated 
monochromatic flux decrement, integrated over the solid angle of 
the source. Triangles shows flux distribution for the run with an 
injection of 1 keV of energy at z = 7. 

ing such requirements, for example the OVRO and BIMA 
interferometers, and they are already being used to detect 
and study SZ clusters. Improvements to these experiments 
should soon reduce the upper bounds shown in Fig. || and 
make the thermal SZ signal we predict detectable in the 
near future. The planned satellite mission Planck Sur- 
veyor with its high sensitivity and multi-frequency capa- 
bility should allow a clean separation of the thermal SZ 
signal from the primary anisotropics, even though its an- 
gular resolution will be limited to less than i ~ 1000. How- 
ever, future interferometric missions at mm- wavelengths, 
like ALMA, will be able to probe much smaller angular 
scales, up to £~10^, and should have enough sensitivity 
to detect both the thermal and kinetic effects. More de- 
tailed discussions of the prospects for extracting the SZ sig- 
nal from upcoming CMB expe riments are given elsewhere 
(Cooray, Hu fc Te gmark |2000| , Holder & Carlstrom |1999 
and Bartlett 200Cl| , among others). 



We have examined the dependence of the SZ signal on 
early energy injection using an "extreme" model where 
0.1 keV to 1 keV of energy was added into all of the gas 
at z = 7. The energy injection heats the gas, but also 
drives it out of small halos, erasing small-scale structure. 
The loss of small scale structures causes a reduction of the 
angular power spectrum at high-^ while the heating of the 
gas enhances the large angle signal (see Fig. 

In the runs without energy injection, most of the mean 
Comptonization is generated at redshifts z < 3, with a 
sm all ta il extending to a redshift of about 5 (da Silva et 
al. 200C). However, in the runs with energy injection, ad- 
ditional Comptonization is provided by hot ionized gas in 
the IGM at still higher redshift. Current FIRAS limits 
on the mean Comptonization already exclude our most 
extreme model that injected 1 keV at z = 7. In future 
work, we will investigate the dependence of the SZ effects 
on additional input physics such as radiative cooling, star 
formation and its associated energy input by supernova 



feedback. 
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